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Abstract. We introduce a new, publicly available Monte Carlo radiative transfer code, Stokes, 
which has been developed to model polarization induced by scattering off free electrons and dust 
grains. It can be used in a wide range of astrophysical applications. Here, we apply it to model 
the polarization produced by the equatorial obscuring and scattering tori assumed to exist in 
£ — , active galactic nuclei (AGNs). We present optical/UV modeling of dusty tori with a curved 
inner shape and for two different dust types: one composition reproduces extinction properties 
of our Galaxy, and the other is derived from composite quasar spectra. The polarization spectra 
enable us to clearly distinguish between the two dust compositions. The Stokes code and its 
i ■ documentation can be freely downloaded from http: / /www.stokes-program.info/j 
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T 7~^ ' 1. Introduction 

' Spectropolarimetric data in the optical/UV range put important constraints on the 
^sO • emission and scattering geometry of AGNs (see e.g. IAntonuccT 2002 1 . For the interpre- 
tation of spectropolarimetric data detailed modeling tools are important. Here, we use 
a new, publicly available radiative transfer code, Stokes (Goosmann & Gaskell 2007), 
to model the obscuring torus of AGNs. The code is based on the Monte-Carlo method 
and allows the simulation of various emission and scattering geometries. Polarization 
due to Thomson and dust (Mie-)scattering is included. Moreover, the code computes 
wavelength-dependent time delays and can thus be used to study polarization reverber- 
i ation (Shoji, Gaskell, & Goosmann 2005). 

^ I 2. Modeling an obscuring torus 

We consider a torus with an elliptical cross-section and centered on a point source. 
The source isotropically emits a flat continuum spectrum between 1600 A and 8000 A. 
The torus half-opening angle is set to do = 30° . The inner and outer radii of the torus 
are fixed at 0.25 pc and 100 pc respectively. The radial optical depth in the equatorial 
plane is 750 for the V-band. The dust models (table [1]) assume a mixture of graphite and 
"astronomical silicate" and a grain radii distribution n(a) oc a" between a m i n and a max . 

The "Galactic dust" model reproduces the interstellar extinction for i?v = 3.1 whilst 
the "AGN dust" parameterization is obtained from quasar extinction curves derived by 
Gaskell et al. 2004. This latter dust type favors larger grain sizes. 



Table 1. Parameterization of the dust models 



Type 


Graphite 


Silicate 








Galactic 
AGN 


62.5% 
85% 


37.5% 
15% 


0.005 fim 
0.005 fim 


0.250 (xxa 
0.200 fj.ui 


-3.5 
-2.05 
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Figure 1. Polarization (top) and flux (bottom) spectra for a centrally-illuminated torus filled 
with "Galactic dust" (left) or "AGN dust" (right). The flux spectra are normalized to the 
illuminating flux and the inclination angle i is measured from the symmetry axis. 

3. Results and discussion 

The angular distributions of the polarization and flux spectra are shown in figure [TJ 
In a face-on view, the central source is seen directly and the obtained polarization is 
low. At obscured inclinations i > 9 , the scattering properties of both dust types lead to 
different results. For "AGN dust", the polarization is lower in the UV but rises quickly 
toward longer wavelengths. Detailed spectropolarimetric observations of dust reflection 
thus effectively constrain the dust composition. For both dust types the polarization 
position angle is oriented perpendicularly to the projected symmetry axis of the object. 

If the torus is filled with "AGN dust" the total flux spectra are wavelength-independent 
above 2500 A at all possible inclinations. For the "Galactic dust" model the scattered 
flux rises gradually toward longer wavelengths. The albedo of both dust compositions 
changes significantly below 2500 A so that less radiation is scattered in the UV. 

Spectropolarimetric data from AGNs contain not only information about the obscuring 
torus, but also include other scattering regions. To obtain a more detailed picture of 
AGN polarization it is therefore necessary to model several, radiatively coupled scattering 
regions self-consistently. While this is beyond the scope of this proceedings note, Stokes 
is capable of solving such problems. 
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